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Abstract. We describe different types of self-trapped optical beams carrying phase 
dislocations, including vortex solitons and ring-like soliton clusters. We demonstrate 
numerically how to create such nonlinear singular beams by the interaction of several 
fundamental optical solitons. Mutual trapping of several solitons can be regarded as 
a synthesis of 'soliton molecules ', and it corresponds to a transfer of an initial orbital 
angular momentum of a system of solitons to a spin momentum of an optical vortex. 
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1. Introduction 

Phase dislocations carried by the wavefront of a light beam are associated with the 
zero-intensity points where the light intensity vanishes. The phase of the wave is twisted 
around such points creating a structure associated with an optical vortex. Optical beams 
with phase dislocations play an important role in linear singular optics jl]. 

In self-focusing nonlinear media, an intense finite-extent laser beam becomes 
localized due to the self-trapping mechanism which can compensate for the beam 
diffraction. The nonlinear self-action of light may result in the formation of stationary 
structures with both intensity and phase remaining unchanged along the propagation 
direction. Such self-trapped stationary structures of light beams are termed spatial 
optical solitons 0. When such solitons have phase singularities, they determine the 
internal structure of the beam; they can be stabilized by the light self-action generating 
nonlinear self-trapped optical beams carrying phase dislocations. Examples of such 
beams include vortex solitons |Sl IH E] with point screw dislocations, multipole vector 
solitons with vr-edge dislocations, and more complicated 'necklace -type beams [11013 
and soliton clusters ^0] with a combination of a screw dislocation at the beam center 
and, generally, -d-edge dislocations, where is the phase jump between neighboring 
peaks in the intensity distribution [TUj . 
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The fundamental optical solitons show a fascinating combination of the properties 
of classical wave-packets together with a number of particle-like properties demonstrated 
in their elastic and inelastic interactions and mutual scattering, when each of the solitons 
preserves its identity. Moreover, the coherent interaction between the solitons depends 
strongly on a relative phase what provides an additional degree of freedom to control the 
interaction. We may draw an analogy between spatial soliton and the 'atom of light', 
and then the soliton collision and interaction can be treated in terms of the effective 
forces acting between these effective 'atoms'. Following this concept, the higher-order 
multi-hump optical beams can be regarded as bound states of 'atoms' trapped by a 
common potential induced in a nonlinear medium. A balance of the interaction forces 
acting between the solitons is the necessary condition for the formation of the soliton 
'clusters' or 'molecules of light'. 

In this paper we investigate the excitation of higher-order beams, including optical 
vortices and soliton clusters, through the inelastic soliton scattering and mutual trapping 
of initially well-separated fundamental solitons, the effect resembling a synthesis of 
'soliton molecules'. In addition, we propose the application of this effect in the context 
of 'soliton algebra' regarding the fundamental spatial solitons as the information 
carriers, and the transformation of an optical pattern induced by the soliton interaction 
as all-optical soliton switching. 



2. Optical vortex solitons and soliton clusters 

Optical vortices were introduced as the first example of a stationary light beam with 
the phase twisted around its core; the twist is proportional to 271 with integer m, the so- 
called topological charge of the phase dislocation ^ . The physical model analyzing the 
evolution of the slowly varying field envelope E is described by the nonlinear Schrodinger 
equation, 

^— + A^E + fiI)E = 0, (1) 

where A_|_ is the transverse Laplacian, and z is the propagation distance measured in 
the units of the diffraction length. Function /(/) describes the nonlinear properties of 
an optical medium, and it is assumed to depend on the total beam intensity, / = \E\ . 
The simplest spatially locahzed solution of equation ((1} carrying a phase dislocation, 
i.e. the vortex soliton, can be written in the form: E{r, 6, z) = A{r)e'^'^^'^^^^ , where A{r) 
and P are the beam amplitude and propagation constant, respectively, while r and 9 
are the polar coordinates in the transverse plane. 

In self-focusing nonlinear media, such vortex beams are the subject of the azimuthal 
modulational instability which result in splitting of the doughnut ring-like structure 
into a certain number of the fundamental solitons. The number of splitters and 
their dynamics are determined by the topological charge of the phase dislocation 
corresponding to the beam angular momentum (see, e.g., |1] and references therein). 
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Figure 1. Intensity and phase distributions for optical vortex soliton and soliton 
cluster composed of four fundamental beams. Note that in terms of the azimuthal 
coordinate 6 = ta.n~^ {y / x) , the vortex phase is given as a linear function mO with 
integer m, while the staircase-like phase of the cluster is a nonlinear phase dislocation. 

The simplest higher-order scalar stationary solutions found for the model (0) are 
a family of the radially symmetric solitons which includes radial modes with nodes in 
form of concentric rings. The important characteristic of these states is the soliton 
spin determined as a ratio of two conserved quantities, the beam angular momentum 
and the beam power. For the vortex soliton, the spin coincides with the topological 
charge m of the phase dislocation carried by the vortex. Because of the condition of 
the field univocally, the topological charge m is quantized and has integer value. The 
fundamental spatial solitons and their higher-order radial states have zero spin. 

Novel types of the higher-order self-trapped optical beams can be introduced as 
azimuthally modulated self-trapped structures in the form of the so-called 'necklace' 
beams [Z1E|- However, it was found that a combination of the edge-phase dislocation 
with TT-out-of-phase neighbor peaks cannot produce a stationary state, and the structure 
becomes slowly expanding T. Such a stabilization is indeed possible for a more 
complicated system including the attraction between several incoherent beams jHl Ej- 
Another approach to this problem is to combine the screw dislocation in the origin 
of a ring-shaped beam with the edge dislocation within the necklace |8]. The screw 
dislocation introduces a centrifugal force to the ring, being also responsible for spiralling 
and mutual repulsion of the solitons in the case of vortex break-up 0, and the edge 
dislocations prevents noise- induced instability break-up of the ring. Because of a nonzero 
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angular momentum, the whole structure rotates with its propagation. As a result, the 
stabilization of the ring-shaped multihump beams requires complex phase distribution 
characterized by fractional value of the soliton spin [Hlini- 

A phase distribution required for the formation of quasi-stationary higher-order 
self-trapped optical beams was found in references where the concept of soliton 
cluster was introduced. In this approach, the azimuthally modulated beam is regarded 
as a bound state of the interacting fundamental solitons. Because of phase-sensitive 
interaction, the requirement of the balance of the interaction forces between the solitons 
determines the beam phase in the form of a staircase-like screw dislocation. Figure 
compares the vortex phase dislocation (left column) with the phase of a four-soliton 
cluster, having well defined 7i/2 steps between the soliton positions (right column). It 
was found ^01 that a radially stable dynamical bound state is formed if these phase- 
jumps satisfy the condition = 2Tim/N , with > 4m being the number of solitons in 
the ring. 

Stability of the soliton clusters has been tested numerically for different nonlinear 
media, including cubic saturating, competing cubic self-focusing and quintic self- 
defocusing, and competing quadratic and cubic self-defocusing nonlinearities (mj. The 
idea has been also extended to higher dimensions, covering the case of the spatio- 
temporal vortex solitons and light bullets. The common outcome of these studies 
is the confirmed robustness of the soliton clusters to random noise and strong radial 
perturbations. In the latter case, the pulsating states viewed as the radial excitations of 
'soliton molecule' have been observed. Nevertheless, soliton clusters are the subject of 
modulational instability, and they are unstable with respect to azimuthal perturbations. 
The remarkable feature of this instability is that the number of the fundamental solitons 
flying off the main ring after the splitting is determined mainly by the topological charge 
m instead of the initial number A^ of solitons, similar to the vortex solitons. For what 
follows, we stress the fact that the conservation of the angular momentum of an optical 
beam in an isotropic medium determines the dynamics of splitters after break-up, so 
that the initial 'spin' angular momentum of vortex or cluster can be viewed as being 
transformed into the orbital momentum of the spiralling splitters |4;^ . 

3. Soliton molecules 

Recent progress in both theoretical and experimental studies of the higher-order optical 
spatial solitons brought the soliton community to the gates of the direct search for 
all-optical soliton-based switching schemes, when the initial data carried by the light 
distribution on the front-face of the nonlinear medium can be processed, in a predictable 
and controllable way, by employing the light self-action effects. One of the examples of 
such an approach is the recently proposed concept of the 'soliton algebra' [TT], based 
on the instability-induced break-up of optical vortices to the controllable number of the 
fundamental solitons. This idea also represents the example of a nontrivial approach 
to the soliton instability, when the symmetry-breaking instability, usually regarded as a 
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Figure 2. Switching from A'' = 6 to A'' = 3 solitons. After reaching the minimum 
radius at 2; = 30, the dynamically unstable three-lobe cluster is formed with its 
subsequent break-up into three fundamental solitons flying away. Note that symmetry- 
breaking instability is accompanied by power and angular momentum loses due to 
nonsolitonic radiation emission. 
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Figure 3. Switching from N = 8 to N = 4 solitons. After reaching the minimum 

radius at z = 18, the ring-like structure creates a single-charged metastable vortex 
which breaks-up subsequently into four fundamental solitons. 



serious disadvantage in using spatial solitons, is employed as a key physical mechanism 
for all-optical soliton switching from a given initial state (optical vortex) to the known 
final state defined by a certain number of fundamental solitons. This approach can be 
generalized to a broad variety of scalar and vector higher-order metastable solitons. 

The symmetry-breaking soliton instability may serve as a physical mechanism for 
all-optical switching with only one disadvantage - it is a one-way process describing the 
transition from an initial higher-order state (a soliton molecule or cluster) to a number 
of simple stable states (dipole-mode and fundamental solitons). Below, we propose 
the opposite process, viewed as the excitation or 'synthesis' of higher-order states from 
predefined number of initially separated solitons, or 'atoms of light', in a nonlinear bulk 
medium. Indeed, introducing molecules of light would not be self-consistent without the 
possibility of mutual trapping of the free atoms or molecule synthesis. 

To demonstrate this phenomenon, we propagate numerically the ring-shaped arrays 
of initially well separated coherently interacting fundamental solitons in a saturable 
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nonlinear medium. Figure |21 shows a characteristic example of a set of six solitons 
which have their relative phases growing in steps of = it/ 3 along the ring, being 
initially directed to collide with each other. This initial condition corresponds to the 
inversion of the instability-induced ring break-up, so that the solitons move towards 
the ring instead of flying away. We observe highly inelastic collision of the solitons 
when they strongly interact overlapping and loosing their identity. Nevertheless, the 
initial phases of the solitons are tilted in such a way that the total phase of the 
beams forms a screw dislocation in the ring origin which prevents a simple fusion of 
all solitons. Instead, the ring-shaped structure is formed. Similar situation is observed 
in figure El with an array of eight solitons and the formation of a metastable vortex 
ring. Due to large amplitude modulations, these intermediate (or metastable) structures 
never form ideal stationary states, and they quickly split off to a set of new isolated 
solitons, with the total number of splitters predefined by the initial conditions. In this 
way, we were able to produce, as a final state, the patterns with different number 
of solitons by changing initial parameters, including the number and phase tilt of 
solitons. Generally, the final number of solitons is determined by the ring instability 
mode with a largest growth rate, and in a saturable medium it is usually twice the 
topological charge j3]. At the same time, we were able to force single-charged 'meta- 
vortex' in figures El and El to split to three and four fundamental solitons, respectively. 
The intermediate meta-state shows complex dynamics of the instability development 
which may continue for several tenths of diffraction lengths. Thus, the whole picture 
of the soliton collision, ring formation, and the ring splitting is somewhat similar 
to the 'delayed-action interaction' recently reported for interacting composite solitons 
carrying nonzero angular momentum ^2j. We note also that, in addition to the known 
transformation of 'spin-to-orbital' angular momentum jlj, in figures |21 and El we observe a 
kind of 'orbital-spin-orbital' transformation. The change of the field momentum shown 
in figure El is about 10% of the initial value, and it occurs only at the break-up stage. 

Highly unstable ring-shaped beams are formed as a result of mutual trapping and 
inelastic interaction of coherent solitons, as shown in figures El and El They represent 
the intermediate steps in the process of the nonlinear switching between the states 
with a different number of solitons. In the context of the soliton algebra and all- 
optical switching, it might be of interest to study in more details and to determine 
the quantitative parameters necessary to obtain the final state with a given number of 
solitons. At the same time, the mutual trapping of coherent solitons is found to be too 
sensitive to the initial perturbations in order to produce the metastable clusters. 

Nevertheless, it is indeed possible to generate higher-order optical beams, i.e. to 
synthesize the soliton molecules, from interacting fundamental coherent and incoherent 
solitons. In figure ^ we show the excitation of a long-lived four-lobe vector cluster. 
We note that the four solitons in the figure jH^b) at z = are directed exactly to the 
center so that there is no azimuthal tilt of their phase. As within the cluster in figure ^ 
each soliton has a phase growing by a jump -d = tt/2 along the ring, and this initial 
phase distribution guarantees the appearance of the screw phase dislocation, even the 
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Figure 4. Excitation of a four-lobe rotating vector cluster, consisting of fundamental 
[row (a) and dotted lines] and modulated [row (b) and dashed lines] mutually incoherent 
components. Solid line: the total angular momentum. Arrows in the row (b) show the 
initial transverse velocities of interacting 'atoms' and rotation of the final 'molecule'. 

corresponding value of the total angular momentum is very small (a solid line in figure^}. 
Because only the total angular momentum is conserved, but not the partial momenta 
of the components, there is a freedom for components to symmetrically exchange the 
angular momentum during the beam propagation. After mutual trapping of all solitons 
at 2; ~ 25 the new-born vector cluster experiences strong radial oscillations between 
the states with maximal (e.g. at 2; = 30 and z = 60) and minimal partial angular 
momenta in the components, see the diagram in figure |3] At the same time, the n/2- 
phase jumps introduced initially between the solitons (edge dislocations) survive these 
strong oscillations, and the whole cluster preserves its structure for a distance exceeding 
100 diffraction lengths. Therefore, we observe the formation of a composite state, the 
soliton molecule, by colliding simple solitons with a nontrivial phase pattern. 

4. Conclusions 

We have studied the scattering and mutual trapping of several fundamental solitons and 
the generation of soliton clusters and vortex solitons - the complex self-trapped states 
of light carrying phase dislocations in the wave front. Inelastic collision of solitons has 
been shown to result in the formation of ring-shaped beams or metastable vortices which 
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subsequently break-up creating different number of fundamental solitons defined by the 
initial conditions. This kind of the soliton delayed-action interaction and the nonlinear 
transformation of the number of fundamental solitons has been analyzed in the context 
of the soliton algebra and all-optical soliton switching. We have also shown that the 
vectorial interaction between the field components provides an additional mechanism of 
the sohton cluster stabihzation. In particular, we have demonstrated the formation of 
vector soliton clusters from colliding solitons, the process which can be regarded as a 
synthesis of 'soliton molecules'. 
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